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INTRODUCTION 

In this analytical work the influence of the Saffman- Taylor instability [1] on flame 
propagation is formulated for computational investigation Specifically, it is of interest to 
examine the influence of this instability as a potential means of eliminating the effect of 
gravitational acceleration on the development of thermoacoustic instability. Earlier experimental 
investigations of thermoacoustic instability [2-6] employed tubes of large circular or annular 
cross-section, such that neither heat loss nor viscosity at the burner walls was of significant 
importance in influencing flame behavior. However, it has been demonstrated recently [7] that 
flames propagating between closely spaced walls, may be subject to long- wavelength wrinkling 
associated with the Saffman-Taylor instability, known to be relevant when a less-viscous fluid 
pushes a more- viscous fluid through a porous medium or between two closely spaced walls. 


FORMULATION 

Conservation equations governing the interaction between a flame and the dynamics of the flow 
through which it propagates have been formulated for conditions under which the activation energy 
characterizing the overall chemical reaction is large and the flame thickness and characteristic 
chemical time are much smaller than the length and time scales characterizing temporal and spatial 
variations of the flow field, respectively [8]. These equations have been solved analytically for freely 
propagating flames with weak curvature and flow-field strain to provide jump conditions for the 
flow field across the flame and an evolution equation governing the motion of the flame [9-12]. 
Recently, the formulation has been extended to account for the presence of acoustic waves in the 
flow that may influence flame dynamics and result in thermo acoustic instabilities [4, 13, 14]. This 
formulation does not take into account, however, influences of heat loss or viscosity at the burner 
walls, either of which could be important if the walls confining flame propagation are sufficiently 
close. In the present work the formulation of Searby and Rochwerger [14] is employed to account 
for the presence of acoustic waves, but with the addition of a term in the momentum conservation 
equation accounting for the effect of viscosity at the burner walls, following the phenomenological 
approach of Joulin and Sivashinsky [7]. The following equations are assumed to govern fluctuations 
of the velocity and pressure field (u, p ), assumed to be small. 
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The independent spatial variables appearing in Eq. (1) are measured in units of the tube diameter, 
annulus width or flat- plate separation distance d , while time is measured in units of d/U L ; where 
U L is the laminar burning speed for the given reactive mixture. Like the independent spatial 
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variables, the dependent flame-surface perturbation a about the location of the mean flame-surface 
plane is also measured in units of d while velocity, pressure and density p are measured in units of 
U L , p U U 2 L and p„, respectively; where p„ is the density of the cold reactant flow. Other 
parameters appearing in the equation are the Prandtl and Richardson numbers defined as Pr =v /a 
and Ri = gd/ul , where v and a are the kinematic viscosity and thermal diffusivity of the mixture, 
respectively; and the ratios e and s of the flame thickness to d and of the thermal conductivity of 
the mixture to its value in the cold reactant flow, respectively. The parameter/ is a constant friction 
factor accounting for the effect of viscosity at the inner walls of the tube or annulus within which, or 
flat plates between which, the flame propagates and is determined by the profile of the mean 
velocity across the flow cross-section. Note that the independent spatial variable x is measured 
normal to the reaction zone of the flame, having been introduced through a transformation from 
laboratory coordinates to a coordinate system that moves with the flame locally [8]. 

The motion of the flame, considered as an interface separating the cold reactants from the hot 
products, is dictated by its advection by the flow and its local rate of propagation into the reactants, 
as determined by thermodynamic and transport properties of the flow and local flow-field strain and 
flame-surface curvature. Small local perturbations about the mean reaction- zone location of a quasi- 
planar, weakly stretched flame (with small flow-field strain and flame-surface curvature) are 
predicted [8-10, 12, 15-18] to obey a linear evolution equation when the overall chemical activation 
energy is assumed to be large; namely, 
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where the velocity fluctuation along the direction of mean flame propagation u and its gradient 
du/dx are evaluated in the reactant flow just upstream of the flame, as denoted. Eqs. (1) and (2) 
comprise a closed set of five equations to be solved for the five flame and flow-field perturbation 
variables (a,p,u ). Since the flow- field equations are valid only in the incompressible regions on 
either side of the flame, jump conditions on the flow variables must be satisfied at the moving flame 
front to ensure continuity of the flow across the flame. The jump conditions are obtained in the 
same manner as that leading to Eq. (2), through integration of the full set of conservation equations, 
including those for energy and reactant species, across the reaction and preheat zones of the flame 
using asymptotic expansions that exploit large chemical activation energies and the disparity between 
the thicknesses of the reactant and preheat zones and between the overall flame thickness and the 
length scales over which order-unity variations of the hydrodynamic reactant and product flow 
occur. 

Evaluation of linear stability proceeds with the assumption of harmonic flow-field and flame- 
surface fluctuations (u_,p,a ) according to e n,+i -i, where k is the wavenumber vector 
characterizing fluctuations in the transverse coordinate plane, measured in units of d ', and Q. 
represents the growth rate (real part) and frequency of fluctuations (imaginary part) measured in 
units of U L Ul . Calculation of the fluctuation amplitudes by solution of Eqs. (1) and (2), subject to 
the jump conditions across the flame, results in a dispersion relation which governs whether the 
fluctuation amplitudes will grow in time or be attenuated by stabilizing influences. The dispersion 
relation, which takes the form 


A(k)a 2 + B(k)a-C(k) = 0, ( 3 ) 

where A , B and C are functions of the magnitude k of the wavenumber vector characterizing 
fluctuations in the transverse coordinate plane and the thermodynamic and transport properties of 
the reactant and product flows, has been derived in earlier studies of the stability of freely 
propagating wrinkled flames with weak curvature and flow- field strain [9-12]. Recently, corrections 
to the functions A , B and C valid through order £ 2 to account for the effect of viscosity at the 
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burner walls has been provided by Joulin and Sivashinsky [7]. However, in this recent work second- 
order influences on flame stability associated with flame stretch were neglected. In the present work, 
the recent corrections have been incorporated into the earlier formulation to provide the following 
forms of the functions A , B and C that account for influences of both viscous damping of the mean 
flow that could give rise to the Saffman-T aylor instability as well as flame-stretch, which may 
provide attenuation of small-scale flame-surface wrinkles. 

A = (g +1)+£&(o -1 ){ Ma — f ^ j ) 
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The new parameters introduced in Eg. (4) are defined as follows, 
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where 0 denotes the ratio of the departure of the temperature from that of the reactants {T -T u ) to 
the temperature difference across the flame (T h -T u ), % is the ratio of the thermal conductivity to 
the specific heat normalized by the same ratio evaluated in the cold reactant mixture, and L<- is the 
Lewis number defined as the ratio of the thermal diffusivity to the mass diffusivity of the limiting 
reactant, where one reactant is considered to be present in excess. The Zel'dovich number (1 is 
defined as the ratio of the activation energy E to the product of the gas constant R and burnt-gas 
temperature T h , normalized by a measure of the temperature difference across the flame, and is 
assumed to be large. The parameter g is simply the ratio of the temperature of the burnt gases to 
that of the cold reactant flow. The dispersion relation as defined in Eg. (3) indicates that flow-field 
and flame-surface perturbations will be unstable and grow exponentially in amplitude when the sign 
of last term of the eguation, C(k ) , is positive. Otherwise, the growth rate Q has no positive real 
part and such perturbations will be attenuated as a result of stabilizing influences. The new terms in 
the dispersion relation accounting for possible influence of the Saffman-T aylor instability appear as 
second-order corrections to coefficients A and B involving the product of the friction factor f and 
the Prandtl number. Note in particular that although the effect of gravitational acceleration acting 
along the mean propagation direction, exhibited through the Richardson number Ri , is stabilizing in 
that the value of C(k) decreases with increasing Ri , the new second-orderviscous term counteracts 
the influence of this stabilizing effect and may cause C(k) to become positive when the tube 
diameter, annulus width or flat-plate separation distance d is sufficiently small. 

The problem of linear stability can be posed alternatively, in terms of the same coefficients A , B 
and C, by considering the dispersion relation as the characteristic eguation associated with a second- 
order evolution eguation governing a linear oscillator. Specifically, the second-order, linear ordinary 
differential eguation 
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has solutions of the form exp(nr) , with the growth rate Q. satisfying the dispersion relation given 
in Eq. (3) resulting from solution of the set of Eqs. (1) and (2) and application of the flow-field jump 
conditions at the flame front. Analysis of Eq. (6) can be fruitful, therefore, in obtaining predictions 
of the stability of flame-surface perturbations once the coefficients A , B and C have been 
determined in the manner described above. Eq. (6) serves also as a point of departure for 
consideration of the influence of acoustic flow-field fluctuations on flame stability, following Searby 
and Rochwerger [14] and Bychkov [131 through the introduction an oscillating acceleration field as 
an additional term in the function C(k). Specifically, replacing Ri in this function with 
Ri-C/ fl cocos(coO, where [J a and go are the amplitude and frequency of axial acoustic velocity 
fluctuations present in the flow measured in units of u, and ujd , respectively, results in the 
following modified form of Eq. (6). 

3 2 oc 

A(k)—r r + B(k)— — [c(&) -D(k)t/ f( oocos(Gor)]oc = 0 

of at , ( 7 ) 

D(k) = ek(a -l)[l-eik(M a — ^j-7)] 

Analysis of Eq. (7) allows the calculation of stability diagrams predicting acoustic-velocity amplitudes 
and wavenumbers for which flame propagation is unstable [13, 14]. The stability diagrams are 
obtained in the same manner as that described in detail by Bychkov [13]. However, predictions made 
with the modified formulation presented herein account for the importance of the Saffman-Taylor 
instability, associated with the disparity between the viscosities of the reactant and product flows, 
when the tube diameter, annulus width or flat-plate separation distance d is sufficiently small. 
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